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Abstract 

Recent measurements of cosmic-ray electron and positron fluxes by PAMELA 
and ATIC experiments may indicate the existence of annihilating dark matter with 
large annihilation cross section. We show that the dark matter annihilation in the 
big-bang nucleosynthesis epoch affects the light element abundances, and it gives 
stringent constraints on such annihilating dark matter scenarios for the case of 
hadronic annihilation. Constraints on leptonically annihilating dark matter models 
are less severer. 



1 Introduction 



Cosmological observations have revealed that about 20 percent of the total energy density 
of the Universe is dominated by the dark matter (DM) [TJ, whose detailed properties 
are still unknown, and many physicists believe that the DM is a kind of stable particle 
appearing in the physics beyond the standard model. Thus, to determine the origin and 
nature of the dark matter in the Universe is one of the most important topics in the 
particle physics, and some methods were proposed for detecting the signals of the DM 
directly or indirectly [2j [3] . One of such methods is that to search for high-energy cosmic- 
rays, including gamma-rays, positrons, anti-protons and neutrinos, which come from the 
DM annihilation in our Galaxy. 

Recent results of the cosmic-ray positron and electron fluxes by the PAMELA satellite 
experiment [I] and the ATIC balloon experiment [5] are now drawing a lot of attention, 
since the steep excess observed by these experiments can be interpreted as an extra 
contribution from the DM annihilation. (For earlier papers, see [6].) However, in order 
to explain these signals the annihilation cross section should be fairly large, as (av) ~ 
(10 —10 -23 ) cm 3 s _1 , which may be achieved by the Sommerfeld enhancement effect [7j. 
This is orders of magnitude larger than the standard value (av) ~ 3 x 10~ 26 cm 3 s _1 , which 
reproduces the observed DM abundance under the thermal freezeout scenario [2]. Then, 
a huge boost factor (Bp ~ 100), which is due to the enhancement of the DM annihilation 
rate due to the clumpy structure of the DM halo, should be introduced in the scenario if 
(av) is not varied in time. However, the DM may be produced nonthermally in other ways, 
such as late-decay of long lived particles [H [9]. Once we give up the thermal freezeout 
scenario, a large annihilation cross section is still allowed. 

The DM with large annihilation rate leads to other observable signatures: gamma- 
rays, anti-protons, and neutrinos!]] In particular, even if the DM only annihilates into 
leptons, the internal bremsstrahlung processes always emit significant amount of gamma- 
rays, and it also predicts a comparable amount of neutrinos. Those may put stringent 
constraints on the DM models [TTJ [T2] . 

Besides these cosmic-ray signals, in this paper we show that the DM annihilation 
also affects the prediction of the big-bang nucleosynthesis (BBN), since it injects high- 
energy particles in the nucleosynthesis epoch, which modify the light elements abundances. 
Such an effect of DM annihilation was pointed out by Jedamzik in Ref. [T3], where the 
modification on the 6 Li abundance due to the hadronic process was emphasized. (See 
also Refs. [H1[T5] for early studies of DM annihilation effects on BBN.) This subject was 
recently studied by four of the present authors in connection with the observed positron 
excess [TB]. In this paper we have performed more systematic studies on the effect of 
DM annihilation on BBN, including the case where the DM annihilates into only leptons, 
motivated by recent results of PAMELA/ATIC. In such a case, the photo-dissociations 
of light elements give constraint. In particular, we found that the 3 He to D ratio gives 
the most stringent constraint on the annihilation cross section, which can be consistent 

1 Indirect detection signatures of non-thermally produced DM were investigated in Ref. [TU] before 
PAMELA and ATIC. 
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with the PAMELA/ATIC results. We also consider the case that the DM annihilates into 
hadrons. Then, the BBN constraint becomes more stringent, and a boost factor larger 
than unity must be introduced in order to account for the PAMELA/ATIC anomalies. 
Therefore, the DM annihilation models as an explanation of the PAMELA/ATIC results 
should be treated carefully so as not to contradict with the BBN constraint presented in 
this paper. 

This paper is organized as follows. In Sec. [2] effects of DM annihilation on BBN 
are explained both for the cases of radiative and hadronic annihilation modes. Current 
status of observations of light element abundances is summarized there. In Sec. [3] the 
resultant constraints on the DM annihilation cross section and also their implications to 
the PAMELA/ATIC observations are presented. Sec. H] is devoted to conclusions and 
discussion. 



2 BBN with Annihilating Dark Matter 
2.1 Light element abundances: theory 

As discussed in the introduction, we consider the BBN scenario with DM particle which 
has sizable annihilation cross section. In such a scenario, the annihilation rate of the DM 
during the BBN epoch may be significant and the abundances of the light elements may 
be affected by energetic particles emitted via the annihilation process. 

For the study of the effects of the DM annihilation during BBN epoch, it is necessary 
to calculate the production rate of energetic particles (i.e., 7, e ± , hadrons, and so on) 
via the annihilation. Denoting the pair annihilation cross section of the DM as (av) , the 
production rate of the energetic particle / is given by 



dt 



-n 



DM 



(av) 



~dE 



(1) 



where fi(E, t) is the energy spectrum of the particle J at a time t, and [dNi / 'dE] ann is the 
energy distribution of / from the single annihilation process. In our study, we assume that 
the pair annihilation is dominantly via s-wave processes, and that (av) is independent of 
time. 

In addition, n^M is the DM number density (at the time t), which is given by 



n DM 



^dmP 



crit 



a(t) 



(2) 



where p cr i t is the present critical density, ttidu is the DM particle mass, and ao and a(t) are 
the scale factors at present and at the time t, respectively. For the density parameter for 
the DM, we take Odm = 0.206 p]. The energy distribution [dNj/dE] ann strongly depends 
on the properties of the DM; it is sensitive to how the DM annihilates. In particular, 
if quarks and/or gluon are produced via the annihilation, we should take account of the 
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Final State 


Eyis/ m D M 


e + e~ 


2.00 




0.70 


r + r - 


0.62 


W + W~ 


0.94 



Table 1: Visible energy carried by e^s and photons after annihilation and subsequent decay, 
£vis, in single DM annihilation process. 



hadronization process in calculating [dNi / dE] ann . We use the PYTHIA package [18] for 
the precise estimation of the distribution functions. 

It should be noted that, as one can see from Eq. ([!]), the production rate of the 
energetic particles (per one DM particle), which is given by ~ nD M (av), is proportional 
to a -3 and is enhanced in the early Universe. This is a significant contrast to the case 
where the DM particle is unstable. The PAMELA and ATIC anomalies can be explained 
if the present production rate of e ± is ~ 10~ 26 s _1 , which may be realized either by the 
annihilation or the decay. If the production of e ± is via the decay, the production rate is 
given by the decay rate of the DM particle and is a constant of time. Then, the production 
rate of the energetic particles is too small during the BBN epoch to affect the abundance 
of light elements. 

Effects of the DM annihilation is classified into (i) photo-dissociation, and (ii) p <-> n 
conversion due to emitted pions and hadro-dissociation. In the following, we summarize 
important points in these effects. 

2.1.1 Photo-dissociation of light elements 

The high energy charged leptons and photons emitted into the cosmic plasma induce 
electromagnetic showers and produce copious energetic photons. These photons destroy 
the light elements ( 4 He, 3 He, D and 7 Li) synthesized in BBN through photo-dissociate 
processes [T9~| [20]. 

The effect of photo-dissociation is determined only by the amount of total visible 
energy E V - 1S of produced particles in the annihilation if E V - 1S ^> 10 MeV [19]. Here "visible 
energy" is defined by the sum of energies carried by e ± s and photons after annihilation and 
subsequent decay. Importantly, E vis depends on the annihilation modes of dark matter. 
We estimate E v - ls using the PYTHIA package; for the final states which will be studied in 
the following, we show the ratio E vm /m-DM in Tabled! (Notice that the ratio E v i S /m-DM is 
independent of mDM once the final state is fixed. In addition, E V - 1S < 2todm because we 
consider pair annihilation processes of the DM.) 

Once E vis is given, the rate of the visible energy injection during the BBN epoch is 
given by 

= ^E vis n 2 DM (av). (3) 
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With the given injection rate, the spectra of the energetic photon and electron induced 
by the DM annihilation are obtained by solving a set of Boltzmann equations, which 
include effects of various radiative processes (photon-photon pair creation, inverse Comp- 
ton scattering, Thomson scattering and so on). Then, we have incorporated the photo- 
dissociation rates in network calculation of BBN and obtained the abundances of the 
light elements. In this calculation we have used the most recent data for nuclear re- 
action rates [2D [221 [221 [2U [22] and estimated the theoretical errors by Monte Carlo 
method [23 [21]. 

The details of the photo-dissociation effects have been discussed, for example, in 
Ref. [20J. The important points of the photo-dissociation effects are summarized as fol- 
lows: 

• T > 10 keV: The emitted high energy photons are quickly thermalized through 
photon-photon process and no significant photo-dissociation of the light elements 
occurs. 

• 1 keV < T < 10 keV: The only D is destroyed by the high energy photons and its 
abundance decreases. The other light elements are not affected. 

• T < 1 keV: All light elements are destroyed by the high energy photons produced 
by DM annihilation. As a result, the abundance of 4 He decreases while D, 3 He and 
6 Li are non-thermally produced. In this case 3 He to D ratio gives the most stringent 
constraint. 

Here we make comments on annihilation onto hadronic particles. Even if DM annihi- 
lates mainly into charged leptons, we expect some emission modes into hadrons via, for 
example, DM + DM — > £ + + £~ + Z* , where denotes charged lepton hereafter. As will be 
seen in next subsection, energetic hadrons, especially nucleons, cause "hadro-dissociation" 
of the light elements and give more serious effect on BBN than high energy electrons and 
photons [23 [23 129], [301 [31] . However, the expected branching ratio into hadrons is ~ 10~ 3 , 
which is too small to give a significant constraint in the present problem. (See Fig. [3j) 

For annihilation into r's, pions are produced from r decays. Approximately 0.9 charged 
pions are emitted by the decay of r lepton in average. As we will see in the next section, 
the effect of pions might be important at around T ~ 1 MeV when the n/p is affected by 
n + 7r + — > p + 7T° and p + 7r _ — > n + 7r°. The p «-> n inter-converting reactions increase 
n/p, which results in larger value of the 4 He abundance. However, we have checked that 
the constraint from the overproduction of 4 He is much weaker than the constraint from 
the photo-dissociation. 

2.1.2 Hadronic effects 

Next we discuss the effects of hadron emissions on the abundances of the light elements. 
As we will see, if the annihilation process has a hadron-emission mode with a net branching 
ratio of > O(10)%, the most stringent constraint is mainly from the hadronic modes. The 
exceptional case is that the emission time is very late (T < 0.3 keV). 
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Once energetic colored particles are emitted by the annihilation, they are hadronized 
and the energetic mesons and baryons are produced. In particular, once energetic nucleons 
are emitted after the hadronization process, they may scatter off the background nuclei 
and induce non-thermal production and destruction processes of light elements. The 
most significant processes are via the scattering off the background 4 He. The effects of 
the non-thermal production of the nuclei A (= D, 3 He, 4 He, 6 Li, and 7 Li) can be taken 
into account by including the following production term into the Boltzmann equations: 

= ^U^nnluiav). (4) 

Here, ^.ann is the number of the produced nucleus A per single annihilation process, 
which depends on the background temperature and the energy spectrum of the emitted 
hadrons. The spectra of the produced hadrons depend on the annihilation mode and the 
mass of the DM. We note here that the mass- dependence is given by ^A,ann oc rn!^ 15 with 
6 < 0.2 [29J. Then, as we see in the following, the constraint from the hadron emission 
approximately scales as oc {av) / m]^. Such processes become effective when the cosmic 
temperature is lower than ~ 0.1 MeV. 

At higher cosmic temperature, another process, which is the p <-> n conversion process 
induced by emitted pions, is more effective. If the cosmic temperature is sufficiently high 
(T > 0.1 MeV), all the emitted hadrons are stopped through electromagnetic interaction 
with background electrons and photons before they scatter off the background nuclei. In 
that case, the stopped hadrons such as charged pions (as well as nn and pp) scatters 
off the background nuclei only through exothermic reactions with their threshold cross 
sections. Then background neutron and proton are inter-converted each other through 
the process like n + n + — > p + ir° andp+7r _ — > n + n . Those processes change the neutron 
to proton ratio even after the freezeout time of neutron. In this case more 4 He and D are 
produced through the non-thermal interconversion processes [HI [32] • Note that unstable 
hadrons such as 7r° decay with short lifetimes (r <C 10~ 8 sec) and disappear well before 
contributing to the inter-converting process. 

We have included the relevant hadro-dissociation processes as well as effects of p <-> n 
conversion into the BBN network calculation and obtained the light element abundances. 
The details are discussed in [251 12"9"1 [50] . and here we just summarize the most important 
features of the hadronic effects: 

• T > 10 MeV: The emitted hadrons cannot change the neutron to proton ratio and 
do not affect any light element abundances at all. 

• 100 keV < T < 10 MeV: 4 He and D are produced through the n <-> p interconversion 
caused by stopped mesons i^) and/or nucleon-antinucleon pairs (pp and nn). 

• 10 keV < T < 100 keV: 4 He is destroyed by energetic nucleons. That produces a 
lot of energetic D, T and 3 He. Among those daughter nuclei, D gives most stringent 
bounds on the amount of the primary nucleons emitted by the annihilation. 
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• 0.3 keV < T < 10 keV: 6 Li is non-thermally produced by the scattering between 
the background 4 He and the energetic T. This gives the severest bound. 

• T < 0.3 keV: The effect from the photo-dissociation can dominate and the 3 He to 
D ratio gives the severest bound when E vis /m DM > 0.1 even if the hadronic mode 
is at 100%. 

2.2 Light element abundances: observations 

Next, we summarize observational constants on primordial abundances of D, 3 He 4 He, 6 Li 
and 7 Li, which we adopt in our study. The errors are presented at 1 a. The subscript "p" 
and "obs" are for the primordial and observational values, respectively. 

The primordial value of D abundance is inferred from observation of high redshift QSO 
absorption systems. We adopt a following observational constraints on the deuterium 
abundance as "Low" value: 

Low (n D /n H ) p = (2.82 ± 0.26) x 10" 5 , (5) 

which is the most-recently reported value derived by taking the weighted mean of six 
observed QSO absorption systems [33] • This value well agrees with the baryon to photon 
ratio suggested by the WMAP 5- year CMB anisotropy observation [I]. However, the six 
measurements have a larger dispersion as expected. In addition, the deuterium is the 
most fragile element. So, in order to derive a conservative constraint, we also adopt the 
highest value among the six measurements as "High" D/H, 

High (n D /n H ) p = (3.98^?) x 10~ 5 . (6) 

To constrain the primordial 3 He abundance, we use an observational 3 He to D ratio 
as an upper bound, which is a monotonically increasing function of the cosmic time. 
(For details, see Refs. [3H El]-) In this study we adopt the newest values of D and 
3 He abundances simultaneously observed in protosolar clouds (PSCs), (n3 He /n H )psc = 
(1.66 ± 0.06) x 10" 5 and (n D /n H ) PSC = (2.00 ± 0.35) x 10~ 5 [35\. Then we get 

(n3 H c/n D )p < 0-83 + 0.27. (7) 

The mass fraction of 4 He is determined by the measurement of recombination lines 
from extragalactic HII regions. In the most recent analysis [36], two values are reported 
by using old and new 4 He-emissivity data, Y p = 0.2472 ±0.0012 and Y p = 0.2516 ±0.0011, 
respectively. Notice that the error presented in [3Sj does not include systematic effects. 
Thus, in this study, we add an error of 0.0040 [3Z] to derive conservative constraint: 

Y p = 0.2516 ± 0.0040. (8) 

(For the systematic uncertainty of the observed 4 He abundance, see also Ref. [58].) 

For 7 Li observed abundance in high temperature metal-poor population II stars is 
considered as primordial. Here, we adopt the most recent value of the 7 Li to hydrogen 



6 



ratio log 10 ( 7 Li/H) obs = —9.90 ± 0.09 given in Ref. [39J. This is close to the value given in 
Ref. [3D] : logi ( 7 Li/H) obs = —9.91 ± 0.10. On the other hand, a slightly larger value has 
been also reported in Ref. [31]: logio( 7 Li/H) b s = —9.63 ± 0.06. The problem is that the 
theoretical value of 7 Li is much larger than the observational one even if we adopt the 
higher value of the observational abundance. The situation becomes worse when we use 
an updated reaction rate of 4 He( 3 He,7) 7 Be [2S1I32I- The observational values in Refs. [39] 
and [31] are smaller than the standard BBN prediction by approximately 0.35 dex and 
0.25 dex, respectively. As for 6 Li abundance, on the other hand, recent observation shows 
that the theoretical value is much smaller than that of the observation, ( 6 Li/ 7 Li) b s = 
0.046 ± 0.022 [33]. These two discrepancies may be collectively called "lithium problem". 
Concerning the inconsistency between the face value of the primordial Li abundance and 
the standard BBN prediction, various solutions have been discussed from the viewpoints 
of both astrophysics [331 133 06] and cosmology [33 [27J [31]. H However, the main purpose 
of this paper is to derive a conservative constraint, so we do not go into the details of 
these models. Instead, assuming some depletion, we add an additional systematic error 
of +0.35 dex into the observational face value n7 Li /nu in our study: 



Notice that this depletion factor D 7 = 0.35 is the systematic error a little larger than that 
allowed from effect of 7 Li depletion by rotational mixing in stars [33] H 

As for a 6 Li constraint, we use (n6 Li /n7 Li ) obs = 0.046 ± 0.022, which was newly- 
observed in a very metal-poor star [33]. We also add a systematic error of +0.106 [50] 
to take into account depletion effects in stars adopting the relation between 7 Li and 6 Li 
depletion factors, D 6 = 2.5D 7 [51], [33] El which leads to A log 10 (77,6 L i/?T.7 Li ) = 0.525 for 
D 7 = 0.35. Then we get a following upper bound: 



3 Constraints on Annihilation Cross Section 

Now, we show the constraints on the annihilation cross section of the DM particle. As 
shown in Eq. ([T]), the amount of energetic particles produced by the annihilation process 
is proportional to the annihilation cross section, and hence the effects on the light element 
abundances become more enhanced as the cross section becomes larger. Since the light 

2 For references, see also [48] for possible astrophysical scenarios to enhance 6 Li abundance, and [49] 
for its observational uncertainties. 

3 Thc depletion factor D 7 is defined as D 7 — A log 10 ( 7 Li/H). The depletion factor for 6 Li is defined 
similarly. 

This relation is obtained for depletion by rotational mixing. If the depletion is caused by atomic 
diffusion [IS], this relation is no longer correct. Instead we should have used D 6 ~ D 7 . However, in this 
case the upper bound on 6 Li/ 7 Li becomes more stringent than given by Eq. (flT)]) . We prefer to use the 
relation for the rotation mixing because we will get more conservative constraint. 



log 10 (m Li /n H ) p = -9.90 ± 0.09 + 0.35. 



(9) 




(10) 
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Figure 1: BBN constraints on the annihilation cross section of DM particles when we take Dj 
= 0.35. The red, cyan, blue, green and brown solid curves represent the upper bounds from 
observational constraints on 3 He/D, D, 6 Li, 7 Li and 4 He, respectively. The name of the element 
is also written by each line. As for D, we plot both upper limits from High and Low values. 
Here we assume that the dark matter annihilates only into e + e~, which means the fraction of 
the visible energy is E v \ s /m = 2. The lines of 4 He, 6 Li and Li do not appear in this figure. For 
reference, the region sandwiched by two dashed lines is allowed by 6 Li when we take D-j = 0. 



element abundances predicted from the standard BBN is consistent with the observations, 
we obtain an upper bound on (av ) in order not to spoil the success of the BBN. 

In this section, we calculate the abundances of light elements for several different DM 
scenarios, taking account of non-thermal production processes discussed in the previous 
section. Then, we compare the theoretical results with observational constraints to derive 
the upper bound on (av). As we will see below, the constraints on the cross section 
depend on the dominant annihilation process; if the DM dominantly annihilates into 
charged leptons and/or photons, the photo-dissociation processes are the most important 
effect, while the hadronic effects becomes more significant when colored particles are 
produced by the annihilation. In the following, we consider several typical cases. 

3.1 Charged-leptonic modes 

The simplest way of explaining the PAMELA and ATIC anomalies is to adopt the possi- 
bility that the DM dominantly annihilates into charged lepton pair: DM + DM — > i + + £~ . 
We start with considering such a possibility. 

We first show the constraints on (av) for the case where DM + DM — > e + + e~ is the 
dominant annihilation process (Fig. [I]). Since E vis oc m DM , the constraints depend on 
(crv)/mT>w The most stringent constraint comes from 3 He/D. 

We can convert the constraints given in Fig. [T] to other cases as far as the effects of 
hadronic decay modes are negligible; since the effects of the photo-dissociation depend 
only on E vis , the upper bounds from the photo-dissociation processes scale as E v ig. In 
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particular, if we use the most stringent upper bound, which is from 3 He/D, we obtain 

(av) < 3.5 x 10- 22 cmV 1 ( J**-) ~* ( ^L) . (11) 
V 1 \2m BM J VI TeV/ 1 ; 

For the case where DM + DM — > fi + + \i~ or r + + r _ is the dominant annihilation 
process, the constraints can be easily obtained from Fig. [1] by using Eq. pip , because the 
photo-dissociation effect on BBN is determined by the total visible energy injection by 
the annihilation process]^ 

For the cases of DM + DM — > fi + + \i~ and r + + r~, the upper bounds are given by 
(av) < 1.0 x 10" 21 cnrV 1 x (m DM /l TeV)" 1 and 1.2 x 10" 21 cmV 1 x (m DM /l TeV) -1 , 
respectively. 

3.2 Hadronic modes 

Next, we consider effects of hadronic modes. Even though we are mainly interested in 
models which can explain the excess of cosmic-ray e ± , it does not necessarily exclude 
models in which significant amount of hadrons are emitted by the annihilation process. 

One of the well-motivated cases is that the annihilation cross section into the W + W~ 
pair is sizable. This is the case when, for example, the Wino in the supersymmetric 
model is the lightest SUSY particle (LSP) and is the DM. As we discuss in the following 
subsection, the PAMELA and ATIC anomalies may be explained in such a case. Then, 
importantly, leptonic decay of the produced W ± is important for the production of cosmic- 
ray e ± , while the BBN constraint is mainly due to the hadronic decay mode. 

In Fig. [2J we show the bound on the annihilation cross section for the case where the 
DM annihilates only into the W + W~ pair. For m^M ^$ 2 TeV, as one can see, the D 
abundance gives the severest constraint on the annihilation cross section if we adopt the 
low value of the observational D; this is due to the hadro-dissociation process of abundant 
4 He. Then, in such a case, we can obtain an approximate constraint on the annihilation 
cross section as a function of the mass: (av) < 10~ 23 cm 3 s _1 (A^/0.8)" 1 (m DM /l TeV) 1 ' 5 , 
with N n being the number of emitted neutrons per single annihilation (~ 0.8 in the case of 
W + W~ emission); here, we have used the approximate relation discussed in the previous 
section: £ A>ann oc m^. 

So far, we have considered cases where energetic charged leptons (in particular, e^) 
are directly produced by the annihilation of the DM or the decay of heavy particle (i.e., 
W^. Such cases are important to explain the behavior of the e ± fluxes observed by the 
PAMELA and ATIC, because e ± produced via the hadronization process is less energetic. 
Consequently, in the light of PAMELA and ATIC, scenarios where the DM annihilates 

5 For T < keV where the photo-dissociation of 4 He is effective, both muon and tau leptons decay well 
before they loose their initial kinetic energies through scattering off the background photons and electrons. 
That is because timescales of the Inverse Compton (IC) and the Coulomb scatterings (CS) are given by 
tic ~ 2.2 x 1CT 10 s (7^±/10 3 )- 1 (T/keV)- 4 (m^±/m e ) 3 and t cs ~ 0.76 x I0 5 s (^±/I0 2 GeV)(T/keV)- 3 
with 7^± being the Lorentz factor, m £ ± the mass, and E^± the energy of charged lepton £ (= /i and r). 
respectively. 
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m (GeV) 

Figure 2: Same as Fig. Q]but for the annihilation into W + W~ pair. Then the fraction of the 
visible energy is E w \ s /m = 0.94. 



only into colored particles (like a quark pair or a gluon pair) are less attractive. However, 
in general, the DM may annihilates dominantly into colored particles. Thus, we also 
comment on constraints on the annihilation cross section in such a case. Here, as an 
example, we consider the case where the DM annihilates dominantly into a bb pair; the 
constraints is shown in the top panel of Fig. [3j The fraction of the visible energy is 
-EVis/moM = 1-04. In this mode, the number of emitted neutrons is an increasing function 
of mDM (0.3-0.5 for 777dm = 100 GeV-1 TeV). For comparison, we also study the case 
that (a7j) D M+DM^b+fe/( crw )DM+DM- + aii = 0.01, with keeping the relation -E vis /777 D M = 1-04; 
the constraints are shown in the bottom panel of Fig. [31 Comparing the top panel with 
the bottom one, we see that the constraints from the photo-dissociation dominates if the 
hadronic mode is less than 1 %. 

3.3 Implication to PAMEL A/ ATIC 

In previous subsections, we have seen that the quantity (av) should be smaller than 
~ O(10~ 24 — 10~ 21 ) cm 3 s _1 in order not to spoil the success of the BBN, if the dominant 
annihilation process of the DM is into e + e~ , r + r~, W + W~, bb, and so on. In fact, 

the upper bound on (av) is of the similar order of that required to explain the PAMELA 
and ATIC anomalies. 

To see this, we plot the positron fraction and the total electron and positron flux 
in Figs. H] and [5], respectively. Here, we adopt the background fluxes of electron and 
positron of Ref. [52] . The annihilation cross section and the DM particle mass (as well 
the propagation model) used in obtaining Figs. H] and [5] is summarized in Table [21 In 
Fig. 1 (Fig. ED, we also show the data points of HEAT [57] and PAMELA @] (BETS 
P], PPB-BETS [53] and ATIC [5])@ From the figures, we can see that DM models with 

6 Recently the HESS collaboration also reported electron flux with E > 600 GeV [SD], which is consis- 
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Figure 3: Same as Fig. Q]but for the annihilation into bb pair (top panel). Then the fraction of 
the visible energy is E v i s /m = 1.04. For comparison, we also plot the case of b + b emission at 
1% in the bottom panel virtually by changing the hadronic branching ratio from 1 to 0.01 by 
hand with keeping the same value of E v i s . 
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the annihilation cross section of (av) ~ 10~ 23 cnr^s" 1 well explains the anomalies in the 
cosmic-ray e ± fluxes. 

To calculate the electron and positron fluxes from the DM annihilation, we have 
adopted the conventional procedure [53]. Approximating the shape of the diffusion zone 
as a cylinder (with half-height L and the length 2R = 40 kpc), we have derived the static 
solution to the following diffusion equation [53] : 

-/(£,£) = K(E)V 2 f(E,x) + —[b{E)f{E,x)\ + Q(E,x), (12) 

where f(E,x) is the electron/positron number density with energy E at the position x, 
which is related to the DM contribution to the electron and positron fluxes as 

*%^(E t S @ ) = ^f{E,2 @ ), (13) 

where c is the speed of light and x Q denotes the position of the solar system. In Eq. (fl2|) . 
K(E) is the diffusion constant, and Q(E, x) is the source term from the DM annihilation. 
For the energy loss rate, we use b(E) = 1 x 10~ 16 (£'/1 GeV) 2 . 

The derived positron and electron fluxes are sensitive to the diffusion constant and 
the half-height of the diffusion cylinder L. The diffusion constant is parametrized as 
K = Kq(E/1 GeV)" 5 . Hereafter we consider the following two propagation models, called 
the MED model: (K , 5, L) = (0.0112 kpc 2 /Myr, 0.70, 4 kpc) and M2 model: (K , 5, L) = 
(0.00595 kpc 2 /Myr, 0.55, 1 kpc), both of which are consistent with observed boron-to- 
carbon ratio in the cosmic-ray flux [55] . 

The source term in Eq. (TT21) is given by 

Q(E,x) = -B F nl Mtnaw (x)(av) 

where riDM,now(^) is the present DM distribution (in particular, in the halo). Dependence 
of the electron and positron fluxes on the DM density profile is very minor. In our 
calculation, we use the isothermal density profile as PDM,now( r ) = m DMKDM,now( r ) = 
0.43(2.8 2 + 8.5 2 )/(2.8 2 + r 2 pc ) GeV cm" 3 , where r kpc is the distance from the Galactic 
center in units of kpc. Furthermore, according to iV-body simulations, the dark matter 
may not be distributed smoothly in our Galaxy and there may be clumpy structures 
somewhere in the Galactic halo. If this is the case, the positron flux may be enhanced 
[56] . In Eq. (I14p such an effect is considered by the boost factor Bp. 

It is found from the figures that the best-fit values of the DM particle mass and 
the annihilation cross section to the PAMELA and ATIC observations depend on the 
annihilation modes. For example, if the DM dominantly annihilates as DM + DM — > 
e + + e~, mpM — 650 GeV and Bp(av) ~ 0.5 x 10~ 23 cm 3 s _1 is preferred to explain the 
anomalieslll As we have seen in the previous subsection, when m DM ~ 650 GeV, the 

tent with ATIC results. 

7 Model parameters used here and hereafter are just for the representative purpose and the fit to 
the observational data is still reasonable even if the parameters are varied up to a few ten's of percent. 
However, the following conclusions are not affected with such changes of parameters. 
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Figure 4: Positron fraction R as a function of positron energy E. (Top) We assume the DM anni- 
hilating into e + e~ with annihilation cross section {av) = 5 x 10~ cm 3 s -1 for ttt-dm = 650 GeV, 
and into /i + yT with {av) = 15 x 10 -24 cm 3 s _1 for wdm = 900 GeV for the MED propaga- 
tion model. (Bottom) We assume DM annihilating into r + r~ with annihilation cross section 
{av) = 4 x 10~ 23 cm 3 s -1 for moM = 1 TeV, and into W + W~ with {av) = 3 x 10~ 23 cm 3 s~ 1 for 
mDM = 800 GeV for the M2 propagation model. Results of PAMELA and HEAT experiments 
are also shown. 
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Figure 5: Total electron and positron flux (times E^) as a function of their energy. Model 
parameters are same as Fig. UJ Results of BETS, PPB-BETS and ATIC are plotted. 
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Final State ttidm (GeV) Bp(av) (cm s ) Propagation Model 



e + e~ 


650 


0.5 x 10" 23 


MED 




900 


1.5 x 10~ 23 


MED 


T + T - 


1000 


4.0 x 10~ 23 


M2 


W + W~ 


800 


3.0 x 10~ 23 


M2 



Table 2: DM particle mass, the product Bp(av), and the propagation model used in Figs. [4] 
and [5j 

BBN constraint requires (av) < 2.3 x 10~ 22 cnrV 1 , and hence the PAMELA and ATIC 
anomalies can be explained even if Bp = 1. For the case where the DM annihilates as 
DM + DM — > W + + W~ , it is remarkable that we need a sizable boost factor (for example, 
> 4), because the observations require (av) = 3 x 10~ 23 cm 3 s for todm = 800 GeV for 
the M2 propagation model, as shown in Figs. 0] and [5j 

4 Conclusions and Discussion 

In this paper we have investigated the effect of the DM annihilation on BBN. Recent ob- 
servations of the cosmic-ray positron/electron excess by the PAMELA/ ATIC experiments 
can be interpreted as the contribution from the DM annihilation. However, this requires 
very large annihilation cross section as (av) ~ 10 -23 cm 3 s _1 unless the boost factor is not 
introduced, and such large annihilation rate can have significant effects on BBN. 

We found that leptonically annihilating DM models are constrained from the upper 
bound on 3 He/D. This limit is compatible with the DM annihilation models for explaining 
the PAMELA and ATIC anomalies for the case of annihilation into e + e~, and t + t~. 

On the other hand, the DM annihilation models with significant amount of hadrons 
are constrained from the observations of D or 6 Li abundances. This gives a stringent 
constraint in some case; the DM models annihilating into IF-bosons for explaining both 
the PAMELA and ATIC anomalies are excluded unless the boost factor larger than unity 
is introduced. 

In this paper we assumed that the DM annihilation cross section is independent of the 
cosmic time. However, our results also have impacts on scenarios with the time-dependent 
cross section. When the DM annihilates through a pole just below the threshold, the an- 
nihilation cross section is enhanced for lower relative velocity of the annihilating DM 
particles |l7j . This enhancement is automatically implemented in the Sommerfeld mech- 
anism [7J [61]. Then, the large boost factor may not be required to realize both the 
thermal freezeout scenario of the DM and the enhanced fluxes consistent with the 
PAMELA/ATIC. However, in such scenarios, it is notable that the annihilation cross 
section during the BBN epoch may be significantly enhanced than that in the present 
Galaxy. Since our results indicate that the DM annihilation cross section in the BBN 
epoch cannot be much larger than that in the Galaxy now, those models are excluded if 
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the enhancement of the cross section during the BBN epoch is sizable. 

Some comments are in order. The DM annihilation in the Galactic center yields 
significant amount of gamma-rays through the cascade decay of the annihilation products 
and/or bremsstrahlung processes, which must be compared with the HESS observation 
[62] . However, the gamma-ray flux significantly depends on the DM density profile, and if 
a moderate profile such as an isothermal one with rather large core radius [53] is chosen, 
the constraint is loosen [llj . 

DM annihilation cross section can also be constrained from the observations of anti- 
protons [Ml [65] in the case of hadronic annihilation mode, and synchrotron radiation 
[XT] [66] produced by the DM annihilation inside the Galaxy. However, these constraints 
more or less suffer from astrophysical uncertainties, such as the density profile of the DM 
halo, the size of the diffusion zone and the distribution of the magnetic field. They lead 
to orders of magnitude uncertainties on the resulting constraints [67J [68]. Taking into 
account these uncertainties, DM annihilation models are consistent with observations 
[651 E]. 

Neutrinos from the DM annihilation at the Galactic center [12] (and possibly from 
the DM trapped inside the Earth [69J) may also be useful as a tool for cross-checking 
signatures of DM annihilation and Super-Kamiokande [70] gives constraints on some DM 
annihilation models. However, it still depends on the DM density profile, though the 
dependence is rather weak. 

The constraints on the DM annihilation cross section from BBN presented in this 
paper is robust, since they do not suffer from large astrophysical uncertainties, compared 
with those DM signatures in the cosmic rays. 
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